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Final  Report 

The  objective  of  the  present  program,  as  stated  in  the  proposal,  was  to  theoretically 
study  the  ignition  and  flame  stabilization  of  hydrogen  in  air  in  high-speed  flows. 
Emphasis  was  on  computational  study  of  the  structure  of  the  chemically  reacting 
high-speed  flows  with  realistic  chemistry,  derivation  of  reduced  mechanisms  of  the 
reaction  system,  and  asymptotic  analysis  of  the  reaction  zone  structure  and  the 
flame  response.  This  objective  was  satisfactorily  accomplished,  as  summarized  in 
the  following  and  reported  in  detail  in  the  four  archival  publications  that  have  so 
far  resulted  from  the  present  program. 

A  numerical  study  of  the  ignition  in  the  supersonic  mixing  layer  with 
detailed  transport  and  chemical  reaction  mechanisms  was  first  conducted.  The 
emphasis  of  the  study  was  on  identifying  the  controlling  chemical  mechanism  in 
effecting  ignition,  on  the  roles  of  thermal  versus  kinetic-induced  ignition  in  which 
heat  release  and  hence  nonlinear  thermal  feedback  are  not  needed  in  initiating 
system  runaway,  and  on  the  consequences  of  imposing  the  conventional  constant 
property  assumptions  in  the  previous  analytical  study.  Results  show  that  the  state 
of  the  hydrogen/oxygen  second  explosion  limit  has  the  dominant  influence  in  the 
system  response  in  that,  for  all  practical  purposes  ignition  is  not  possible  when  the 
air  stream  temperature  is  lower  than  the  crossover  temperature,  even  allowing  for 
viscous  heating.  On  the  other  hand,  when  the  air  stream  temperature  is  higher 
than  the  crossover  temperature,  the  predicted  ignition  distance  indicates  that 
ignition  is  feasible  within  practical  supersonic  combustion  engines.  For  the  latter 
situations,  the  ignition  event  is  initiated  by  radical  proliferation  and  hence  rimaway 
instead  of  thermal  runaway.  Furthermore,  it  was  demonstrated  that  the  use  of  the 
similarity  Blasius  velocity  at  the  leading  edge  can  generate  a  high  temperature  spot 
which  is  convected  downstream  and  artificially  facilitates  ignition.  This  work  was 
reported  in  Publication  No.  1,  which  was  appended  in  the  Annual  Report  for  1995- 
96. 

The  second  study  was  an  asymptotic  analysis  of  the  ignition  of  the 
hydrogen/air  supersonic  mixing  layer  using  reduced  reaction  mechanisms.  Two 
distinct  reduced  mechanisms  for  the  high-temperature  and  low-temperature 
regimes  were  used  depending  on  the  characteristic  temperature  of  the  reaction  zone 
relative  to  the  crossover  temperature  at  which  the  reaction  rates  of  the  H+O2 
branching  and  termination  steps  are  equal.  Each  regime  also  involves  two  distinct 
analyses  for  the  hot-stream  and  the  viscous-heating  cases,  depending  on  the  relative 
dominance  of  external  and  internal  ignition  energy  sources.  These  four  cases  were 
been  analyzed  separately,  and  it  was  shown  that  the  present  analysis  successfully 
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describes  the  ignition  process  by  exhibiting  turning  point  and  thermal  runaway 
behavior  in  the  low-temperature  regime,  and  radical  branching  followed  by  thermal 
runaway  in  the  high-temperature  regime.  Results  for  the  predicted  ignition 
distances  were  then  mapped  out  over  the  entire  range  of  parameters,  showing 
consistent  behavior  with  the  previous  one-step  model  analysis.  Furthermore,  it  was 
demonstrated  that  ignition  in  the  low-temperature  regime  is  controlled  by  a  large 
activation  energy  process,  so  that  the  ignition  distance  is  more  sensitive  to  its 
characteristic  temperature  than  that  in  the  high-temperature  regime.  The  ignition 
distance  was  also  found  to  vary  non-monotonically  with  the  system  pressure  in  the 
manner  of  the  well-known  hydro  gen /oxygen  explosion  limits,  thereby  further 
substantiating  the  importance  of  chemical  chain  mechanisms  in  this  class  of 
chemically-reacting  boundary  layer  flows.  This  work  was  reported  in  Publication 
No.  2,  which  was  appended  in  the  Annual  Report  for  1995-96. 

Further  improvements  in  the  reduced  mechanisms  were  conducted,  and 
ignition  in  the  counterflow  has  been  studied.  The  first  study  was  concerned  with 
ignition  in  the  high  pressure  situations  of  the  third  limit.  In  the  analysis, 
simplifications  to  the  chemical  mechanism  were  formally  derived,  and  criteria 
determining  the  accuracy  of  these  simplifications  were  identified.  It  was 
demonstrated  that  the  chemical  steady-state  approximations  of  O  and  OH  are 
accurate  for  a  wide  range  of  conditions,  while  the  steady-state  approximations  of  H 
and  HO2  were  both  found  to  be  inaccurate  in  certain  regimes.  Using  these 
approximations,  a  reduced  mechanism  was  derived  and  used  to  explain  the 
ignition  behavior  of  H2  and  02-  It  was  shown  that,  for  both  the  premixed  and 
nonpremixed  configurations,  ignition  occurs  through  the  effect  of  heat  release  on 
the  Arrhenius  term  of  reaction  H2O2  +  M  ->  20H  +  M,  and  expressions  were  derived 
which  were  able  to  qualitatively  predict  the  relationship  between  the  pressure, 
temperature,  and  strain  rate  at  ignition.  The  difference  in  ignition  temperature 
between  the  premixed  and  nonpremixed  configurations  is  associated  with  the  effect 
of  reactions  between  H  and  HO2.  For  the  premixed  case,  the  rate  of  destruction  of  H 
through  H+O2+M  ->  HO2+M  is  much  slower  due  to  the  lower  concentration  of  O2 
in  the  reactions  zone,  thus  reactions  between  H  and  HO  2  play  a  greater  role  in  the 
ignition  process.  This  work  was  reported  in  Publication  No.  3,  and  attached  herein 
as  Appendix  A. 

A  second  analysis  was  performed  for  the  situations  of  high  temperature  and 
low  pressure.  Starting  with  seven  elementary  reactions,  a  reduced  mechanism  was 
derived  assuming  steady  state  for  the  O  and  OH  radicals.  An  algebraic  ignition 
criterion  was  derived  using  this  mechanism.  This  criterion  successfully  explains  the 


behavior  analogous  to  the  first  and  second  explosion  limits.  A  bifurcation  analysis 
was  then  performed  to  clarify  the  ignition  behavior.  This  analysis  demonstrated 
that  an  ignition  turning  point  can  occur  solely  through  the  interaction  of  radical 
species  with  no  contribution  from  heat  release.  The  source  of  this  turning  was 
found  to  be  the  reaction  H  +  HO 2 ->  20 H,  confirming  results  from  numerical 
calculations.  The  regimes  in  which  abrupt  or  monotonic  transition  to  an  ignited 
state  were  re-calculated  including  the  effect  of  this  reaction.  This  work  was  reported 
in  Publication  No.  4,  which  is  attached  herein  as  Appendix  B. 

Publications 

1.  "A  numerical  study  of  ignition  in  the  supersonic  hydrogen/air  laminar 
mixing  layer,"  by  M.  Nishioka  and  C.K.  Law,  Combustion  and  Flame,  Vol.  108,  pp. 
199-219  (1997). 

2.  "Ignition  in  the  supersonic  hydrogen/air  mixing  layer  with  reduced 
mechanisms,"  by  H.G.  Im,  B.T.  Helenbrook,  S.R.  Lee  and  C.K.  Law,  Journal  of  Fluid 
Mechanics,  Vol.  322,  pp.  275-296  (1996). 

3.  "Ignition  of  hydrogen  and  oxygen  in  counterflow  at  high  pressures,"  by  B.T. 
Helenbrook  and  C.  K.  Law,  Twenty-Sixth  Symposium  (International)  on 
Combustion,  The  Combustion  Institute,  Pittsburgh,  PA.  pp.  815-822  (1997). 
[Appendix  A]. 

4.  "Theory  of  thermal-  and  radical-induced  ignition  of  counterflowing  hydrogen 
versus  oxygen  at  high  temperatures,"  by  B.T.  Helenbrook,  H.G.  Im  and  C.K.  Law, 
Combustion  and  Flame,  in  press.  [Appendix  B]. 
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IGNITION  OF  HYDROGEN  AND  OXYGEN  IN  COUNTERFLOW 

AT  HIGH  PRESSURES 

B.  T.  HELENBROOK  and  C.  K.  LAW 
Department  of  Mechanical  and  Aerospace  Engineering 
Princeton  University,  Princeton,  N]  08544 


Ignition  of  hydrogen  and  oxygen  in  the  “third  limit"  was  studied  in  counterflow  for  both  premixed  and 
nonpremixed  systems  with  activation  enei^  asymptotics.  The  results  of  the  analysis  were  then  compared 
to  numerical  calculations  for  further  verification  of  the  findings.  In  the  analysis,  simplifications  to  the 
chemical  mechanism  were  formally  derived,  and  criteria  determining  the  accuracy  of  these  simplifications 
were  identified.  It  was  demonstrated  that  the  chemical  steady-state  approximations  of  O  and  OH  are 
accurate  for  a  wide  range  of  conditions,  while  the  steady-state  approximations  of  H  and  HO2  were  both 
found  to  be  inaccurate  in  certain  regimes.  Using  these  approximations,  a  reduced  mechanism  was  derived 
and  used  to  explain  the  ignition  behavior  of  H2  and  Og.  It  is  shown  that,  for  both  the  premixed  and 
nonpremixed  configurations,  ignition  occurs  through  the  effect  of  heat  release  on  the  Arrhenius  term  of 
reaction  15,  H2O2  4-  M  —>  20H  +  M,  and  expressions  were  derived  that  are  able  to  qualitatively  predict 
the  relationship  between  the  pressure,  temperature,  and  strain  rate  at  ignition.  The  difference  in  ignition 
temperature  between  the  pre mixed  and  nonpremixed  configurations  is  associated  with  the  effect  of  reac¬ 
tions  between  H  and  H02-  For  the  premixed  case,  the  rate  of  destruction  of  H  through  H  4-  O2  4-  M 
HO2  4r  M  is  much  slower  due  to  the  lower  concentration  of  O2  in  the  reaction  zone,  thus  reactions  between 
H  and  HO2  play  a  greater  role  in  the  ignition  process. 


Introduction 

Hydrogen-oxygen  combustion  is  of  interest  be¬ 
cause  it  plays  an  underlying  role  in  the  combustion 
of  all  hydrocarbon  fuels  and  because  hydrogen  is  a 
practical  fuel.  This  interest  has  motivated  a  substan¬ 
tial  amount  of  research  on  its  oxidation  and  com¬ 
bustion  characteristics.  Indeed,  work  done  by  Lewis 
and  von  Elbe  [1],  Baldwin  and  Mayer  [2],  Baulch  et 
al.  [3],  Wamatz  [4],  and  Yetter  et  al.  [5]  has  resulted 
in  a  reliable  description  of  the  chemical  kinetics  of 
hydrogen  oxidation,  which  has  enabled  further  stud¬ 
ies  of  the  behavior  of  hydrogen -oxygen  mixtures  in 
various  physical  environments  [6-8].  The  counter¬ 
flow  mixing  configuration  has  been  extensively 
adopted  for  experimental,  numerical,  and  analytical 
studies  of  hydrogen -oxygen  combustion  [9-14]  be¬ 
cause  it  embodies  many  of  the  physical  processes 
seen  in  practical  applications.  Results  from  such 
studies  have  been  able  to  assess  the  interaction  of 
chemistry  with  convection  and  diffusion. 

Numerical  and  experimental  studies  of  hydrogen- 
oxygen  ignition  in  the  counterflow  configuration 
have  mapped  the  response  of  the  system  over  a  wide 
range  of  conditions  for  both  premixed  and  nonpre¬ 
mixed  systems.  Similar  to  the  homogeneous  explo¬ 
sion  limits  of  hydrogen-oxygen  mixtures,  three  dis¬ 
tinct  ignition  limits  in  such  heterogeneous, 
mixing-affected  systems  were  found  [11].  In  the  first 
limit,  ignition  is  governed  by  a  balance  between 


branching  and  diffusive  loss  of  H,  while  the  second 
limit  is  determined  by  a  balance  between  branching 
and  destruction  of  H  through  H  +  O2  +  M  — »  HO2 
+  M.  In  the  third  limit,  the  branching  rate  of  the 
hydrogen  radical  is  not  sufficiently  fast  to  cause  ig¬ 
nition,  and  ignition  proceeds  through  a  path  involv¬ 
ing  hydrogen  peroxide.  In  addition,  numerical  sim¬ 
ulations  have  led  to  the  suggestion  of  various  skeletal 
mechanisms  and  simplified  global  mechanisms  ap¬ 
plicable  to  each  limit  [9,11]. 

For  nonpremixed  ignition,  asymptotic  analyses 
have  been  performed  in  the  first  and  second  limits 
[12-14]  as  well  as  the  third  limit  [12].  The  analysis 
of  Ref.  12  for  the  third  limit,  however,  was  based  on 
a  reduced  mechanism  obtained  from  empirical  ob¬ 
servations  of  numerical  results.  As  such,  only  limited 
insight  was  obtained.  Therefore,  the  primary  objec¬ 
tive  of  the  present  investigation  is  to  provide  a  rig¬ 
orous  and  comprehensive  analysis  of  the  third  igni¬ 
tion  limit,  for  both  premixed  and  nonpremixed 
systems,  and  through  this  analysis  gain  fundamental 
understanding  of  the  ignition  process.  In  this  anal¬ 
ysis,  we  will  derive  a  global  mechanism  from  a  re¬ 
duced  set  of  reactions  that  has  been  shown  to  ac¬ 
curately  predict  ignition  temperatures  by  previous 
numerical  work  [9,11]  and  has  also  been  extensively 
tested  in  the  course  of  the  present  investigation  for 
both  premixed  and  nonpremixed  systems.  By  deriv¬ 
ing  the  global  mechanism  asymptotically,  we  can  find 
criteria  governing  the  accuracy  of  this  mechanism, 
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and  thus  avoid  some  of  the  inadequacies  that  can  be 
associated  with  determining  the  mechanism  strictly 
from  numerical  results.  The  formulation  and  various 
aspects  of  the  analysis  are  sequentially  presented  in 
the  following. 


Formulation 

Ignition  of  H2-O2  will  be  studied  in  the  counter¬ 
flow  configuration.  We  are  especially  interested  in 
finding  a  relation  between  the  strain  rate  and  the 
temperature  and  pressure  of  the  flow  at  the  point  of 
ignition.  Two  distinct  cases  will  be  considered.  In  the 
first,  the  nonpremixed  case,  a  cold  jet  of  inert-di¬ 
luted  H2  impinges  on  a  hot  jet  of  air.  In  the  second, 
the  premixed  case,  a  cold  jet  of  premixed  H2  and  O2 
impinges  on  a  hot  inert  jet.  Due  to  similarities  in  the 
ignition  behavior  between  these  two  cases,  the  anal¬ 
ysis  for  both  can  be  considered  in  parallel. 

We  will  use  numerical  results  to  guide  the  asymp¬ 
totic  analysis  and  to  verify  its  findings.  These  results 
have  been  calculated  using  a  Newton  iteration  tech¬ 
nique  developed  by  Kee  et  al.  [15]  and  modified  for 
the  counterflow  configuration  by  Nishioka  et  al.  [16]. 
This  code  solves  a  quasi-one-dimensional  problem 
based  on  the  similarity  approximation.  The  calcula¬ 
tions  were  done  using  mixture-averaged  diffusion 
coefficients,  potential  flow  boundary  conditions,  and 
the  reaction  mechanism  given  by  Ref  5. 

For  simplifications,  we  shall  assume  that  the  flow 
is  potential  described  by  /w  ==  (-pjaxa,  pjir), 
where  x  is  the  coordinate  normal  to  the  stagnation 
plane,  p^  the  density  of  the  hot  stream,  a  the  strain 
rate  of  the  flow,  and  a  ==  1,  2  for  planar  and  axisym- 
metric  flows,  respectively.  Furthermore,  Cp,  X,  and 
pDi  for  all  species  were  taken  to  be  constants,  while 
the  Lewis  numbers  for  all  species  were  assumed  to 
be  unity  except  for  H2  and  H,  for  which  the  values 
of  0.3  and  0.2  were  used.  With  the  above  approxi¬ 
mations,  the  governing  equations  in  the  nondimen- 
sional  spatial  coordinate,  x/JX/Cpp^,  become: 


^  (lx  dx^ 


pai(J 


in 


(1) 


^  1 

^  dx  Lei 


W. 

pjia 


(2) 

where  T  =  c„fWJQ,  Y,  =  t.WyW,.  and  q,„  = 
-  v"J}i°/Q. 

For  nonpremixed  ignition,  the  boundary  condi¬ 
tions  were  set  to  describe  air  flowing  against  a  mix¬ 
ture  of  H2  and  No- 

T  =  t1  =  yo2.^ 


y.  =  0  (i  #  O2,  N2);  ^ 


T  =T^-  Yh2  = 

Y,  =  H2.N2);  x-^-00  (3) 


where  P  -  -  T_«.  Furthermore,  we  shall  use 

=  0.1  and  f  =  300  K  in  the  calculations. 
For  pre mixed  ignition,  a  mixture  of  oxygen  and  hy¬ 
drogen  flows  against  a  stream  of  nitrogen,  and  only 
the  boundary  conditions  for  oxygen  and  nitrogen 
were  changed  so  that  for  x  — >  00,  yo2  ”  ^  for  .t 
—00,  Yq^  =  Y'o2,-<»- 

The  frozen  solutions  Tf  and  y,'y  are  given  in  the 
form  Cl  +  C2edc{xjLei/2),  where  Ci  and  C2  are  cho¬ 
sen  to  satisfy  the  boundary  conditions  given  above. 
This  completes  formulation  of  the  problem.  An  anal¬ 
ysis  of  the  above  equations  and  comparison  of  their 
predictions  with  the  numerical  results  follows. 


Large  Activation  Energy  Approximations 

The  first  step  in  the  analysis  is  to  identify  the  char¬ 
acteristic  chemical  rate  at  ignition.  We  will  show  that 
for  both  premixed  and  nonpremixed  systems,  the 
only  consistent  choice  for  this  rate  is  that  of  reaction 
15.  Since  this  reaction  has  an  activation  energy  that 
is  much  larger  than  the  ambient  thermal  energy  of 
the  system,  we  expect  the  analysis  to  follow  proce¬ 
dures  developed  for  one-step  large  activation  energy 
ignition  problems  [17]. 

Inserting  T  =  TV  +  in  the  Arrhenius  term  of 
reaction  15  reveals  that  the  proper  scale  for  e  is 
RT^E^  15  and  that  the  solution  domain  is  divided 
into  three  asymptotic  zones.  Near  the  hot  boundary, 
T„  ~  Tf  is  less  than  0(£),  and  the  Arrhenius  term  of 
reaction  15  is  nearly  constant.  In  the  “reaction”  zone, 
the  fuel  concentration  is  much  larger  than  near  the 
hot  boundary,  while  the  Arrhenius  term  decreases 
exponentially  due  to  an  0(£)  decrease  in  the  frozen 
temperature,  and  in  the  downstream  zone,  reaction 
15  is  essentially  frozen. 

Since  the  reaction  zone  is  located  where  T*  -  7y 
is  of  0(e),  we  can  expand  this  relation  to  reveal  that 
the  physical  coordinate  in  this  zone  can  be  expressed 
as  X  =  Xr  +  <^/Xr,  whcre  changes  in  the  Arrhenius 
term  with  are  0(1)  and  x^  is  a  moderately  large 
number.  Thus,  reaction  will  occur  near  the  hot 
stream  boundary  in  a  narrow  zone  of  width  lAv 
around  the  location  x^.  In  Fig.  1,  we  have  plotted  the 
H2O2  production  rate,  the  H2O2  mass  fraction,  and 
the  temperature  at  ignition  from  the  numerical  re¬ 
sults.  Examining  the  figure  reveals  that  H2O2  is  pro¬ 
duced  in  a  relatively  narrow  zone  near  the  hot 
stream  boundary  from  which  it  diffuses  and  convects 
outward.  In  this  figure,  we  have  also  shown  the  po¬ 
sition  and  width  of  the  reaction  zone  found  from  the 
asymptotic  expression.  This  estimation  agrees  rather 
well  with  the  numerical  results  which  confirms  that, 
typical  of  many  combustion  problems,  the  location 
of  the  reaction  zone  is  determined  by  the  sensitive 
nature  of  the  reaction  rate  to  temperature.  This  is 
true  for  both  the  premixed  and  nonpremixed  prob¬ 
lems. 
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Fig.  1.  Ignition  profiles  at  a  pressure  of  5  atm  and  a 
strain  rate  of  10  s“  ^ 


Applying  the  above  approximations,  we  can  re¬ 
write  the  governing  equations  in  this  zone  as  follows: 


h  Y  Y 

jc  '  I  Ko  ^  m.2 

ci^  E  ))i 

^  .  JL^  ^ 

di  Ui  r 


where 


K.  = 


I 

exp(E^., ,,0/E,,  is) 


(Tax; 

L(Tf) 

liisiTJ ' 


and 


(4) 


(5) 


(6) 


These  equations  determine  the  ignition  behavior 
for  both  the  premixed  and  nonpremixed  cases.  In 


the  above  equations,  A  is  the  Damkohler  number 
based  on  the  convective-diffusive  rate  in  the  ignition 
zone  and  the  rate  of  reaction  15.  Since  this  reaction 
determines  the  characteristic  chemical  rate  of  igni¬ 
tion,  we  expect  the  Damkohler  number  to  be  0(1) 
at  ignition.  Analysis  of  these  equations  will  confirm 
this  expectation,  as  well  as  determine  the  Damkohler 
number  s  dependence  on  the  remaining  parameters 
in  the  system. 


Reaction  Zone  Chemistry 

The  next  step  in  the  analysis  is  to  examine  the 
chemistry  occurring  in  the  reaction  zone.  Examining 
the  reactions  in  Table  1  reveals  that  several  species 
are  consumed  at  a  rate  that  is  much  faster  than  the 
convective-diffusive  rate  in  the  ignition  zone.  The 
radicals  O,  OH,  and  H  all  fall  into  this  category.  Since 
consumption  of  these  radicals  occurs  quickly,  any 
radical  produced  will  be  consumed  before  diffusive 
loss  can  play  a  role  in  this  process.  Making  the  chem¬ 
ical  steady-state  approximation,  we  ignore  the  con¬ 
vective-diffusive  terms  and  balance  the  chemical 
terms  which  gives  an  algebraic  expression  for  the 
radical  cxjncentration. 

For  the  O  radical,  the  consumption  reaction  is  re¬ 
action  2  which  has  a  rate  of  Ak2  Yh,/  relative  to  the 
diffusive  rate  of  the  system.  Sinchez  et  al.  [14]  have 
shown  that  for  ignition  in  the  first  limit,  the  steady- 
state  approximation  of  O  is  inaccurate.  In  the  third 
limit,  however,  the  characteristic  ignition  rate  is 
much  slower,  and  chemical  steady-state  of  the  O  rad¬ 
ical  is  an  accurate  approximation.  The  mass  fraction 
of  O,  therefore,  can  be  found  by  balancing  the  chem¬ 
ical  production  and  destruction  rates  of  O.  Similarly 
for  the  OH  radical,  the  rate  of  consumption  through 
reaction  3  relative  to  the  diffusive  rate,  is 

1  such  that  this  radical  is  also  in  steady  state. 
Because  the  mass  fraction  in  the  reaction  zone  is 


TABLE  1 

Skeletal  mechanism 


No. 

Reaction 

A 

n 

E„ 

R1 

H  +  OH  +  O 

1.92E  14 

0.00 

16.44 

R2 

O  -h  Ho OH  H 

5.08E  +  04 

2.67 

6.29 

R3 

OH  +  Hi  ^  HiO  +  H 

2.16E  +  08 

1.51 

3.43 

R9 

H  +  Oi  +  M  ^  HOi  +  M 

6.17E  +  19 

-1.42 

0.00 

RIO 

HOi  +  H  -*  Hi  +  Oi 

6.63E  13 

0.00 

2.13 

RlOb 

Hi  +  Oi-*  HOi  +  H 

1.93E  +  14 

0.00 

59.61 

Rll 

HOi  +  H  ZOH 

1.69E  4-  14 

0.00 

0.87 

R14 

HOi  HOi  — ^  HiOi  O2 

3.02E  +  12 

0.00 

1.39 

R15 

HjO,  +  M  ^  20H  +  M 

1.20E  +  17 

0.00 

45.50 

RlTb 

HOi  +  Wj  ->  HiOi  +  H 

2.81E  +  11 

0.54 

24.09 

Units  are  mol,  cm*'^,  s,  K,  and  kcal;  k  -  AT^exp  {-EJRT) 
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the  same  for  both  the  premixed  and  nonpremixed 
fuel  streams,  these  approximations  hold  for  both 
configurations. 

Using  these  approximations  in  the  reaction  mech¬ 
anism,  we  find  the  well  known  H  branching  path 
through  reactions  1,  2,  and  3,  which  occurs  at  a  rate 
of  2kiY^YQ^f.  In  the  third  limit,  this  rate  is  slower 
than  the  consumption  rate  of  H  through  reaction  9, 
k^Y^Yo^jlCi^],  Therefore,  the  combined  result  of 
these  reactions  is  consumption  of  ff  at  a  rate  of 
2^1  (/c  ~  1)702/ relative  diffusive  rate  of  the 

H  radical,  where  k  is  the  ratio  of  the  rate  of  reaction 
9  to  the  branching  rate,  kg[Ci^]/2ki,  which  depends 
on  the  pressure  through  [Ca/],  Since  ki  is  of  0(10^ 

10^),  we  expect  the  chemical  steady-state  approx¬ 
imation  of  H  to  be  accurate  when  k  is  larger  than 
one  and  there  is  sufficient  in  the  reaction  zone. 
For  the  nonpremixed  case,  the  O2  mass  fraction  is 
0(1)  in  the  reaction  zone,  and  this  approximation  is 
accurate  even  for  pressures  at  which  k  is  nearly  one. 
In  the  premixed  configuration,  the  reaction  zone  is 
in  the  hot  inert,  and  oxygen  must  diffuse  into  the 
reaction  zone  from  the  fuel  stream.  For  /c  -  1  of 
0(1),  the  O2  concentration  in  the  reaction  zone  is 
enough  to  make  the  chemical  steady-state  approxi¬ 
mation  of  H  fairly  accurate,  but  because  the  oxygen 
concentration  rapidly  decreases  upstream  of  the  re¬ 
action  zone,  convective-diffusive  loss  of  the  H  radical 
may  still  have  some  influence.  For  larger  a:,  this  in¬ 
fluence  is  reduced,  and  the  approximation  becomes 
increasingly  accurate. 

Plotting  the  local  destruction  rate  of  these  radicals 
along  with  their  convective-diffusive  rate,  as  in  Ref 
11,  verifies  the  accuracy  of  these  approximations. 
Figure  2  was  generated  from  the  numerical  results 
for  a  premixed  fuel  stream  at  a  pressure  of  7.4  atm, 
a  strain  rate  of  13  s"^  and  at  a  temperature  of  1020 
K  which  is  within  1  K  of  the  ignition  point.  At  these 
conditions,  k  is  approximately  2.5.  The  figure  shows 
that  for  the  H  radical,  even  for  these  near-worst-case 
conditions,  the  destruction  rate  is  nearly  ten  times 
larger  than  the  convective-diffusive  rate,  thus  chem¬ 
ical  production  is  mainly  balanced  by  destruction. 
For  the  O  and  OH  radices,  the  destruction  rates  are 
more  than  an  order  of  magnitude  larger  than  the 
convective-diffusive  rate,  making  the  chemical 
steady-state  approximation  very  accurate.  This  veri¬ 
fies  that  these  approximations  do  predict  the 
behavior  of  the  numerical  system. 

We  can  use  the  preceding  approximations  to  de¬ 
rive  the  following  reduced  mechanism  which  is  valid 
only  for  ignition  in  the  third  limit: 


Step  I 


{HO2  +  H2O) 


Step  III 


K  + 

1/2 

K  - 

1 

1 

K  — 

K  + 

1/2 

K  - 

1 

2a:  - 

'  1 

K  - 

{H2  +  O2)  H2O2  + 


„  2k  -  1/2  ^ 
H2  ^  ^ 


K  -  I 


HO2  + 


K  -  1 


HoO 


Step  IV  k:(4  -  2a)  -  (1  -h  a) 

2(7-1) 

k{2  -  4fl)  +  1  +  a 


Ho  + 


2{k  -  1) 
2Ka  -  a  —  1 


1 


HOo 


O2  4- 
2k  -  a 


HoO 


2{k  -  1)  - 


R14  HO2  +  HO2  H2O2  +  O2 

Steps  I,  II,  and  III  occur  at  the  rates  of  reactions  15, 
17b,  and  10b,  respectively.  Due  to  the  branching  of 
the  H  radical  and  consumption  through  reaction  9, 
each  H  or  OH  radical  produced  by  these  reactions 
results  in  the  production  of  k/(k  -  1)  HO2  radicals. 
Step  IV  is  the  combined  effect  of  reactions  1 1  and 
10  where  a  is  the  ratio  of  the  rates  of  these  reactions, 
which  is  approximately  0.2.  This  reaction  is  relatively 
slow,  and  will  be  considered  only  as  a  correction  to 
the  leading  order  solution.  Steps  I,  II,  III,  and  re¬ 
action  14  determine  the  primary  ignition  mechanism 
of  the  system. 


Ignition  Mechanism 

The  remaining  analysis  focuses  on  the  system  evo¬ 
lution  as  the  Damkohler  number  is  increased  to  the 
turning  point.  To  understand  this  evolution,  we  will 
begin  by  examining  the  behavior  of  the  HOo  radical. 
In  the  above  reduced  mechanism,  HO2  has  two  pos¬ 
sible  branching  paths.  The  most  apparent  is  through 
step  II,  which  has  a  rate  proportional  to  HOo  and 
results  in  the  production  of  1/{k  -  1)  HO 2  radicals. 
The  second  path  is  through  a  coupling  of  steps  I  and 
II.  For  pressures  at  which  k  is  moderately  larger 
than  one,  this  path  is  the  faster  route,  and  corre¬ 
sponds  to  the  following  pathway: 


HO2  4-  Ho  — ^  H2O2  4*  H 


reaction  17b 


H2O2  +  M  20 H  reaction  15 

20H  +  2H2  — >  2H  +  2H2O  reaction  3 

37/  4-  3O2  4-  3M  3//O2  4-  3M  reaction  9 


(8) 
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Fig.  2.  Destruction  rates  and  convective-difFusive  rates 
of  O.  OH,  and  H  at  7 A  atm,  13  s~^  and  1020  K,  which  is 
within  one  degree  of  the  ignition  temperature. 


I/a[s] 

Fig.  3.  Maximum  HOz  mass  fraction  versus  l/a  for  var¬ 
ious  initiation  rates  and  a  pressure  and  temperature  of  5 
atm  and  970  K. 


When  the  Damkohler  number  is  of 

O  iJ(K  -  the  branching  rate  of 

this  path  will  be  the  same  order  as  the  convective- 
diffusive  rate.  For  Damkohler  numbers  less  than  this 
value,  the  HOo  concentration  is  determined  mainly 
by  a  balance  between  the  initiation  reaction,  RlOb, 
and  the  convective-diffusive  terms.  For  Damkohler 
numbers  larger  than  this  value,  the  rate  of  the 
branching  reaction  becomes  faster  than  the  convec¬ 
tive-diffusive  rate,  and  the  HO2  concentration  in¬ 
creases  to  a  level  at  which  the  dominant  balance  is 
between  the  production  of  HO2  through  this  branch¬ 
ing  path  and  destruction  of  HO2  by  reaction  14.  To 
verify  this,  the  maximum  mass  fraction  of  HOo  has 
been  plotted  versus  l/a  for  several  initiation  rates  in 
Fig.  3.  The  calculations  were  performed  at  a  pres¬ 


sure  of  5  atm  and  a  temperature  of  970  K.  At  these 
conditions,  the  asymptotic  expression  predicts  that 
at  l/a  of  order  0.004  s,  the  transition  between  these 
behaviors  will  occur.  On  the  figure,  the  two  regimes 
are  identifiable.  For  l/a  less  than  approximately  .007 
s,  the  sensitivity  to  initiation  is  apparent  while  above 
this  value,  the  maximum  mass  fraction,  which  typi¬ 
fies  the  variation  of  the  profile  in  the  reaction  zone, 
becomes  fairly  insensitive  to  the  different  initiation 
rates.  This  confirms  that  the  branching  path  is  cru¬ 
cial  in  determining  the  HO2  concentration  at  the  ig¬ 
nition  point. 

For  Damkohler  numbers  much  larger  than  the 
value  at  which  the  branching  increase  of  HO2  takes 
place,  we  can  neglect  the  convective-diffusive  terms 
relative  to  the  branching  rate.  This  leads  to  the 
chemical  steady-state  approximation,  2u;i4  = 
[2kwi^  +  uJiTb  +  {2k  —  f)-  Since 

the  production  rate  of  HO2  through  step  I  exceeds 
that  from  steps  II  and  III,  we  can  ignore  the  second 
two  terms  on  the  right-hand  side  of  this  equation 
and  solve  for  the  concentration  of  HO2: 


^H02  ~ 


ki^iK  -  1) 


^//202 


(9) 


For  this  expression  to  be  correct,  the  Damkohler 
number  at  ignition  must  be  much  greater  than  that 
at  which  HO2  branching  first  occurs. 

Our  final  goal  is  to  determine  the  mechanism  that 
causes  the  ignition  turning  point.  Using  the  previous 
findings  we  can  rev^rite  the  equations  for  H2O2  and 


<t>: 


^f/202  ,  ^^^H202  -  _A; 

•  7  wO  “I 


a:  -  1/2 


■ 


15JW2O2 


K  -  1 


kl7b^H2/ 


kisK 


k^iK  -  1) 


'  W2O2 


(10) 


d(}>  <P(f> 


kisK 


ki4(K  -  1) 


^  H202 


(11) 


Examining  the  above  equations  reveals  that  there 
are  two  sources  for  H2O2  and  cf).  The  first  is  a  branch¬ 
ing  path  with  a  rate  proportional  to  the  concentra¬ 
tion  of  H2O2.  This  path  is  simply  step  I  followed  by 
the  conversion  of  HO2  to  H2O2  through  reaction  14, 
and  will  cause  a  branching  ignition  when  A/(k  —  1 ) 
is  0(1).  The  second  path  corresponds  to  step  II  with 
HO2  converted  to  H2O2  by  reaction  14.  Because  the 
rate  of  this  path  depends  on  the  mass  fraction  of 
H2O2  to  the  one-half  power,  this  path  will  only  cause 
a  quadratic  increase  in  the  concentration  of  H2O2  as 
the  Damkohler  number  is  increased.  A  turning  point 
can  occur  only  when  the  effect  of  heat  release  on 
the  reaction  rates  is  considered.  Determining  the 
amount  of  H2O2  produced  by  this  path  and  inserting 
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Fig.  4.  Ignition  pressure  versus  temperature  for  various 
density-weighted  strain  rates. 


it  into  equation  11  results  in  the  magnitude  of  the 
temperature  perturbation.  Both  paths  are  important 
in  determining  the  ignition  Damkohler  number,  al¬ 
though  the  heat  release  from  path  I  is  larger.  Using 
this  path  to  determine  an  approximate  scaling  for  the 
ignition  Damkohler  number,  we  arrive  at  the  follow¬ 
ing  expression; 


ekuiK  -  IP  1^^ 

qiK{K  -  l/2)2(kl76'^H2.r)^- 


(12) 


where  the  reaction  rate  constants  are  evaluated  at 
T*.  and  Y^2.r  Is  the  mass  fraction  in  the  reaction  zone. 
Comparing  this  expression  to  the  criterion  for  igni¬ 
tion  through  branching  reveals  that  ignition  through 
heat  release  occurs  well  before  branching  of  H2O2 
becomes  important  for  all  k;  therefore,  heat  release 
is  the  dominant  ignition  path.  For  /c  —  1  of  0(1), 
this  expression  is  an  0(1)  quantity  which  confirms 
that  reaction  15  determines  the  characteristic  chem¬ 
ical  rate  at  ignition.  This  is  in  agreement  with  the 
results  found  by  Lee  and  Law  [12]  for  k  much  larger 
than  one.  However,  we  have  now  shown  that  for  k 
of  0(1),  heat  release  is  still  dominant  over  the 
branching  path  in  determining  the  ignition  Dam¬ 
kohler  number. 

Using  the  leading-order  scaling,  we  can  also  de¬ 
termine  a  correction  of  the  form  id  =  -d®  4-  .  .  . 

caused  by  the  effects  of  step  IV.  Assuming  that  the 
rate  of  reaction  1 1  is  small,  we  arrive  at  the  following 
correction  to  the  heat  release  terms  on  the  right- 
hand  side  of  equation  11: 


{-2Ka  +  1  H-  g) 


(1  -h  d«)2(2K:dO  +  dO  4-  1) 

4-  qiv 


A«(2/cd«  4-  /c)  . 

(>-  -  1  -IS)  * 


(13) 


from  which  we  can  determine  the  influence  of  step 
rV  on  the  ignition  Damkohler  number.  The  first  two 
terms  in  this  expression  result  from  the  effect  of  step 
IV  on  the  rates  of  steps  I  and  II.  The  last  term  is  the 
additional  heat  release  derived  from  this  step.  Eval¬ 
uating  this  expression,  we  have  found  that  for  most 
K,  the  effect  of  reaction  11  on  the  rates  of  steps  I 
and  II  is  dominant  over  the  heat  release  from  this 
step.  Near  k  =  3,  the  effect  on  the  reaction  rates 
becomes  very  small  because  step  IV  neither  pro¬ 
duces  nor  destroys  HO^  at  this  k.  For  k<3,  because 
of  the  branching  of  the  H  radical,  step  IV  actually 
begins  to  produce  additional  HO2.  The  ma^itude 
of  this  effect  depends  on  the  rate  of  reaction  11. 
Because  the  H  radical  concentration  is  inversely  pro¬ 
portional  to  its  consumption  rate,  2ki(/c  “•  l)Yo2,f, 
step  IV  will  be  much  faster  near  /c  ==  1  than  for  large 
K.  In  addition,  since  this  rate  is  slower  in  the  pre¬ 
mixed  configuration  due  to  the  lower  oxygen  con¬ 
centration.  step  IV  will  have  a  greater  influence  on 
the  premixed  results.  Therefore,  this  expression  will 
explain  some  of  the  differences  between  the  pre¬ 
mixed  and  nonpremixed  numerical  results. 

The  final  step  in  the  analysis  is  to  check  that  the 
steady-state  approximation  of  HO2  made  previously 
is  a  valid  approximation.  To  do  this,  we  have  calcu¬ 
lated  the  ratio  of  the  Damkohler  number  given  by 
equation  12  to  the  Damkohler  number  found  for 
HO 2  branching  for  various  conditions.  As  previously 
mentioned,  if  this  ratio  is  »  1,  the  branching  rate  of 
HO2  is  much  faster  than  its  convective-diffusive  rate, 
and  the  chemical  steady-state  approximation  is  ac¬ 
curate.  The  results  of  these  calculations  reveal  that 
this  ratio  is  only  between  2  and  10  for  most  condi¬ 
tions  and  that  it  decreases  with  increasing  ignition 
temperature.  Therefore,  the  chemical  steady-state 
approximation  is  not  very  accurate,  and  convective- 
diffusive  loss  of  the  HO2  radical  is  important  in  de¬ 
termining  the  ignition  temperature,  especially  at 
high  temperatures.  However,  as  a  first  approxima¬ 
tion,  chemical  steady-state  of  HO2  still  predicts  the 
correct  qualitative  behavior. 


Comparison  to  Numerical  Results 

In  Fig.  4  we  have  plotted  the  ignition  pressures 
and  temperatures  for  various  density-weighted  strain 
rates  calculated  from  the  numerical  results,  along 
with  the  behavior  predicted  by  equations  10  and  1 1, 
including  the  correction  from  expression  13.  Here, 
instead  of  holding  the  strain  rate,  a,  fixed,  we  hold 
pja  fixed,  as  in  Ref  11.  is  defined  as  the  ratio  of 
pji  to  the  value  at  ac  =  1  for  1  atm.  Holding  d  con¬ 
stant  is  more  appropriate  because  it  fixes  the  scale 
of  the  problem,  and  thus  isolates  the  chemical  ef¬ 
fects.  On  the  plot,  numerical  calculations  for  both 
the  premixed  and  nonpremixed  configurations  are 
shown.  Examining  the  nonpremixed  ignition  results 
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reveals  that  the  scaling  from  the  asymptotic  analysis 
predicts  the  ignition  hmits  fairly  well.  At  high  pres¬ 
sures,  the  ignition  hmits  are  well  predicted  by  bal¬ 
ancing  the  chemical  time  of  reaction  15  with  the 
convective-diffusive  time  in  the  reaction  zone  as 
shown  by  equation  12,  For  k  approaching  1,  the  ig¬ 
nition  temperature  decreases  due  to  the  increased 
effect  of  the  hydrogen  branching  path.  Below  k  = 
1,  the  H  radical  is  no  longer  in  steady  state,  and  the 
system  transitions  to  the  first  limit  behavior. 

Examining  the  premixed  ignition  results,  we  see 
that  the  leading-order  solution  along  with  the  cor¬ 
rection  from  expression  13  does  qualitatively  predict 
the  correct  behavior  The  difference  in  ignition  tem¬ 
perature  between  the  premLxed  and  nonpremixed 
results  can  be  understood  by  examining  this  correc¬ 
tion.  For  large  /c,  the  difference  in  ignition  temper¬ 
ature  between  the  premixed  and  nonpremixed  re¬ 
sults  at  constant  k  is  small,  about  10  K.  At  these 
conditions,  step  IV  is  relatively  slow  because  the 
large  consumption  rate  of  H  causes  its  concentration 
to  be  low.  As  K  decreases,  the  difference  in  ignition 
temperature  between  the  premixed  and  nonpre¬ 
mixed  streams  increases,  reaching  a  maximum  of 
nearly  100  K,  This  is  again  due  to  step  IV.  For  de¬ 
creasing  /c,  the  H  concentration  increases,  increasing 
the  rate  of  this  step  and  thus  making  ignition  more 
difficult.  For  /c  <  3,  the  ignition  temperature  of  the 
premixed  stream  decreases  with  decreasing  k,  rap¬ 
idly  approaching  the  ignition  temperature  of  the 
nonpremixed  stream.  As  previously  mentioned,  for 
K  ^3,  step  IV  no  longer  consumes  HO2,  but  begins 
to  produce  HO 2-  This  along  with  the  0(1)  effect  of 
H  branching  on  the  ignition  temperature  results  in 
the  rapid  decrease  in  ignition  temperature  as  k  ap¬ 
proaches  one. 

The  premixed  case  shown  in  Fig.  4  also  shows  a 
relatively  large  discrepancy  between  the  analytical 
and  computational  results  near  k  =  1.  At  these 
states,  step  IV  becomes  an  0(1)  effect.  Not  only  does 
this  affect  the  leading-order  analysis,  but  it  also  af¬ 
fects  the  steady-state  approximation  of  H.  Since  re¬ 
action  11  is  a  branching  path  for  H,  branching 
through  this  path  can  exceed  the  consumption  rate 
causing  convective-diffusive  loss  of  H  to  become  im¬ 
portant.  For  K  of  0(1),  this  effect  results  in  the  dif¬ 
ference  in  ignition  temperatures  between  the  anal¬ 
ysis  and  the  numerical  results,  and  as  k  approaches 
1,  causes  the  transition  to  first  limit  behavior. 


Conclusions 

Ignition  of  counterflowing  hydrogen  and  oxygen 
was  studied  asymptotically  for  ignition  in  the  third 
limit,  with  the  goal  of  understanding  the  ignition 
mechanism  of  this  system.  Because  of  the  temper¬ 
ature-sensitive  nature  of  reaction  15,  we  found  that 
ignition  occurs  in  a  localized  kernel  near  the  high- 


temperature  boundary.  In  addition,  the  analysis  re¬ 
vealed  that  the  radicals  H,  OH,  and  O  are  in  chem¬ 
ical  steady  state  because  each  is  consumed  at  a  rate 
that  is  much  faster  than  the  diffusive  rate  of  the  sys¬ 
tem.  Using  steady-state  approximations  for  these 
radicals,  a  reduced  mechanism  was  derived  which 
simplified  the  remaining  analysis.  In  analyzing  the 
kinetics  of  HO2,  it  was  found  that  HO2  has  a  branch¬ 
ing  path  that  when  balanced  with  reaction  14  deter¬ 
mines  the  concentration  of //O2  at  the  ignition  point. 
Equations  governing  H2O2  and  the  temperature  per¬ 
turbation,  (f>,  were  derived  which  showed  that  in  the 
third  limit,  ignition  occurs  through  the  effect  of  heat 
release  on  reaction  15,  and  were  able  to  predict  the 
dependence  of  the  ignition  temperature  on  strain 
rate  and  pressure  for  both  nonpremixed  and  pre¬ 
mixed  fuel  streams.  The  difference  in  ignition  tem¬ 
perature  between  these  two  configurations  was  as¬ 
sociated  with  the  effect  of  reactions  of  H  with  HOo. 
For  the  premixed  configuration,  the  H  concentration 
is  much  larger,  such  that  reactions  of  H  with  other 
radical  species  has  a  greater  influence  on  the  ignition 
temperature. 
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Abstract 

Ignition  of  hydrogen  and  oxygen  in  counterflow  was  studied  using 
asymptotic  methods  for  temperatures  above  the  “crossover  tempera¬ 
ture”  .  Starting  from  eight  elementary  reaction  steps,  a  reduced  mech¬ 
anism  was  formally  derived  which  revealed  that  ignition  in  the  high 
temperature  regime  is  characterized  by  a  balance  between  branch¬ 
ing,  termination  and  transport  of  the  hydrogen  radical  within  the 
reaction  zone.  Additionally,  an  algebraic  ignition  criterion  was  de¬ 
rived  which  predicts  the  ignition  state  as  a  function  of  the  parameters 
defining  the  system.  This  criterion  successfully  explains  the  behav¬ 
ior  analogous  to  the  “first”  and  “second”  explosion  limits  observed  in 
homogeneous  hydrogen-oxygen  mixtures.  A  bifurcation  analysis  was 
then  performed  to  clarify  the  ignition  behavior.  It  was  found  that 
there  are  three  distinct  sub-regimes  of  the  high-temperature  ignition 
regime.  Specifically,  at  ignition  pressures  for  which  the  rate  of  the 
termination  reaction,  H -I-O2  + Af  — >  H02-\- M ,  is  comparable  to  the 
branching  rate,  abrupt  ignition  occurs  through  the  nonlinear  effect  of 
the  reaction  H  +  HO2  ^OE  on  the  net  branching  rate.  At  mod¬ 
erately  low  pressures,  production  of  HO2  becomes  negligible,  and 
the  HO2  chemistry  no  longer  controls  the  ignition  process.  At  these 
pressures,  abrupt  ignition  is  caused  by  the  effect  of  heat  release  on 
the  branching  rate.  Finally,  at  very  low  pressures  depletion  of  oxygen 
in  the  reaction  zone  prevents  an  abrupt  ignition,  and  a  monotonic 
transition  to  an  ignited  state  occurs.  The  behavior  predicted  by  the 
asymptotic  analysis  is  in  agreement  with  predictions  from  numerical 
calculations  with  full  chemistry.  In  addition,  the  asymptotic  analysis 
isolated  the  different  physical  mechanisms  underlying  the  ignition  be¬ 
havior  and  determined  parameters  governing  the  transition  between 
these  mechanisms. 


‘Present  address:  Center  for  Turbulence  Research,  Stanford  University,  Stanford,  CA 
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Introduction 


Because  hydrogen-oxygen  combustion  is  of  interest  for  its  application  in  su¬ 
personic  propulsion,  for  its  fundamental  role  in  hydrocarbon  combustion, 
and  for  its  potential  as  a  clean  burning  energy  source,  there  has  been  sub¬ 
stantial  research  done  on  its  oxidation  and  combustion  characteristics.  An 
extensive  amount  of  work  (Baulch  ct  ul.^  1972j  ^Varnatz,  1984,  Yetter  ct  cil.j 
1991)  has  been  performed  on  determining  the  elementary  chemical  mech¬ 
anism  which  governs  hydrogen  oxidation,  and  has  produced  a  reliable  de¬ 
scription  of  the  associated  oxidation  kinetics.  This  chemical  mechanism  has 
since  been  used  to  study  the  behavior  of  hydrogen-oxygen  combustion  nu¬ 
merically  in  various  physical  configurations.  Specifically,  numerical  investi¬ 
gations  have  been  performed  on  homogeneous  mixtures  (Maas  and  Warnatz, 
1988)  and  then  extended  to  more  complicated  situations  involving  diffusive 
transport  (Balakrishnan  et  o/.,  1995;  Darabiha  and  Candel,  1992;  Kreutz  et 
ai,  1994;  Vlachos  et  al.,  1993).  These  investigations  have  shown  that  this 
chemical  mechanism,  although  simple  compared  to  most  hydrocarbons,  can 
produce  surprisingly  complex  behaviors. 

All  of  the  previous  studies  have  reported  that  an  important  parameter 
in  hydrogen-oxygen  ignition  is  the  “crossover  temperature”,  T*.  This  is  the 
temperature  at  which  the  hydrogen  branching  rate  is  equal  to  the  recom¬ 
bination  rate  of  the  hydrogen  radical  with  oxygen  at  a  given  pressure.  For 
temperatures  above  T*,  the  branching  rate  exceeds  the  termination  rate, 
and  ignition  is  governed  by  the  fate  of  the  hydrogen  radical.  For  temper¬ 
atures  below  r*,  the  branching  rate  of  the  hydrogen  radical  is  not  suffi¬ 
ciently  fast  to  cause  ignition,  and  ignition  takes  place  through  an  alternate 
path  involving  hydrogen  peroxide.  The  crossover  temperature  therefore  di¬ 
vides  hydrogen-oxygen  ignition  into  two  distinct  regimes  each  having  its 
own  unique  behavior.  This  study  will  focus  on  ignition  which  occurs  at 
temperatures  greater  than  the  crossover  temperature. 
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Numerical  and  asymptotic  studies  of  the  high  temperature  regime  have 
shown  that  the  ignition  behavior  can  vary  depending  on  parameters  such 
as  the  temperature,  pressure,  reactant  concentration,  and  the  physical  ge¬ 
ometry.  For  example,  a  numerical  investigation  of  the  counterflow  config¬ 
uration  by  Balakrishnan  et  al.  (1995)  has  shown  that  for  ignition  temper¬ 
atures  much  greater  than  T*,  abrupt  ignition  does  not  occur.  Instead,  a 
monotonic  transition  to  a  strongly  burning  state  is  observed.  This  unique 
behavior  of  high-temperature  ignition  has  been  further  investigated  in  a 
bifurcation  analysis  by  Sanchez  et  al.  (1994).  They  have  confirmed  that 
when  the  temperature  is  sufficiently  large,  a  supercritical  bifurcation  is  ob¬ 
served  as  the  Damkohler  number  increases,  implying  a  monotonic  transition 
to  the  strongly  burning  state.  Furthermore,  the  analysis  has  also  revealed 
that,  as  the  temperature  approaches  the  crossover  value,  thermal  feedback 
changes  the  behavior  from  a  supercritical  to  a  subcritical  bifurcation,  and 
thus  the  conventional  turning-point  ignition  is  resumed.  This  study  there¬ 
fore  suggests  a  more  refined  classification  of  the  ignition  behavior  for  high 
temperatures,  depending  on  whether  the  temperature  is  high  enough  to 
suppress  the  effect  of  thermal  feedback. 

A  numerical  study  of  the  same  problem  by  Kreutz  et  al.  (1994)  has  found 
that  when  the  temperature  of  the  hot  air  stream  is  moderately  higher  than 
T*,  the  ignition  turning-point  is  attained  purely  by  a  nonlinear  coupling 
between  radicals,  without  consideration  of  heat  release.  Specifically,  this 
coupling  is  effected  by  an  additional  branching  step  H  +  HO2  —*■  20 H, 
and  constitutes  an  additional  sub-regime  in  the  high-temperature  ignition 
chemistry. 

The  present  study  attempts  to  resolve  all  of  the  possible  factors  that 
are  responsible  for  these  distinguishing  ignition  characteristics  observed  in 
high-temperature  ignition.  In  doing  so,  we  first  notice  that  the  previous 
asymptotic  studies  have  either  used  reduced  mechanisms  which  have  as¬ 
sumed  the  steady-state  approximation  for  the  hydroperoxyl  radical  (Sanchez 
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et  al.  1994)  or  omitted  the  elementary  step  H  +  HO2  — ^  20 H  from  the 
reduced  mechanism  (Lee  and  Law,  1994),  thereby  ruling  out  the  possibility 
of  radical  coupling  observed  by  Kreutz  et  al.  (1994).  The  present  analy¬ 
sis,  therefore,  starts  with  the  derivation  of  a  reduced  mechanism  from  eight 
elementary  reaction  steps  based  on  asymptotics.  Using  this  mechanism, 
we  shall  perform  a  bifurcation  analysis  analogous  to  that  of  Sanchez  et  al. 
(1994).  By  exploiting  the  asymptotic  limit  of  small  Lewis  number  for  the 
hydrogen  radical,  we  obtain  a  simple  algebraic  expression  for  the  ignition  cri¬ 
terion  which  elucidates  the  physical  interpretation  of  the  high-temperature 
ignition  phenomena  in  terms  of  all  system  parameters.  It  is  found  that  high- 
temperature  hydrogen-oxygen  ignition  exhibits  three  distinct  sub-regimes, 
namely  radical-radical  coupling,  radical-thermal  feedback,  and  monotonic 
transition,  as  the  temperature  increases  from  the  crossover  temperature. 
We  shall  also  identify  parameters  governing  the  transitions  between  these 
three  cases. 


Formulation 


The  problem  considered  is  a  counterflow  configuration  in  which  a  stream 
of  cold  hydrogen  impinges  on  a  hot  stream  of  oxygen-nitrogen  mixture.  A 
potential  flow  solution  is  used  to  describe  the  flow  field  in  which  the  mass 
flux  is  given  by  pv  =  ( — p^axcr^  poar).  In  this  expression,  x  is  the  similarity 
coordinate  normal  to  the  stagnation  plane,  po  is  the  density  of  the  incoming 
oxidizer  stream,  a  is  the  strain  rate  of  the  flow,  and  c  is  a  constant  which 
is  equal  to  one  for  planar  flow  and  two  for  axisymmetric  flow.  We  further 
assume  constant  and  equal  values  of  Cp  for  all  species,  and  constant  values 
of  A  and  pDi.  Unity  Lewis  number  for  all  species  is  cilso  assumed  except  for 
H2  and  H,  for  which  the  value  of  0.2  was  used.  Using  the  nondimensional 
spatial  coordinate,  xj ,JDh la,  the  governing  equations  become; 


dT  ^  dPT 

X—  + 

dx  dx^ 


m 


(1) 
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(2) 


^  LencPYj 
^  dx  ^  Lei 


Wo, 


UC 


where  Vq^  is  the  mass  fraction  of  oxygen  in  the  incoming  oxidizer  stream, 
T  =  c^fWoJQYS,,  Yi  =  YiWoJYS,Wi,  =  E.(<n.  "  ^JKIQ,  and 

=  ni[poYm,iYSjWo,]AmT^expi-Ta,m/T). 


Equations  (1)  and  (2)  are  subject  to  the  following  boundary  conditions: 


T  =  Too,  Yo,  =  1,  Yi  =  0(i^  O2,  N2)]  X  00, 

T  =  Too-/3,  Yh,=Y^^,  Yi  =  0iij^H2,N2);  x^-oo, 


where  ^  =  Too  —  T_oo  with  f-00  maintained  at  300K  for  this  problem.  When 
chemical  reaction  is  frozen,  the  temperature,  hydrogen,  and  oxygen  mass 
fractions  are  described  by: 

T/  =  Too  —  ^erfc(x/y2i/eif),  Yh,^  =  -^erfc(x/\/2),  and 

2  ^  (4) 

Yo,,f  =  1  -  -erfc(x/^2i/eH), 

which  will  be  used  as  the  basic  solution  for  the  asymptotic  analysis  to  be 
presented  in  the  next  section.  We  remark  that  the  boundary  conditions  are 
applied  at  finite  distances  away  from  the  stagnation  plane  and  are  taken 
to  the  limit  of  large  distances.  This  is  to  avoid  mathematical  difficulties 
arising  from  imposing  the  boundary  conditions  at  infinity. 

To  proceed  with  the  analysis,  we  need  to  define  the  chemical  mechanism 
of  hydrogen  oxidation.  This  mechanism  is  adequately  described  by  consid¬ 
ering  eight  species  and  nineteen  reversible  reactions  (Yetter  et  ai,  1991). 
Numerical  results  (Trevino,  1991;  Kreutz  and  Law,  1993)  have  shown  that 
at  temperatures  above  T*,  the  eight  elementary  reaction  steps  shown  in  ta¬ 
ble  1  adequately  determine  the  ignition  behavior.  Hydrogen  peroxide  does 
not  play  a  role  in  this  mechanism  because  at  temperatures  greater  than  T* , 
H2O2  reactions  are  too  slow  to  influence  the  ignition  process.  Once  the 
problem  is  solved  one  can  easily  verify  that  the  neglected  reactions  have 
negligible  rates. 
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Asymptotic  Analysis  and  Results 


Anticipating  that  the  ignition  process  is  governed  by  the  fate  of  the  hydrogen 
radical,  we  normalize  the  reaction  rate  constants  in  the  problem  by  the 
rate  constant  of  the  branching  reaction,  Rl,  at  Too-  Since  at  ignition  only 
small  perturbations  to  the  frozen  solutions  are  expected,  we  expand  the 
temperature  as  T  =  Tf  +  (p,  where  <f)  Tf.  This  results  in  governing 
equations  of  the  form; 


where 


dtp 

X—  +  Len 

ax 


<P<P 

dx^ 


m 


^  Lend^Yj 
^  dx  Lei  dx"^ 


(5) 

(6) 


A  =  ^  ^iexp(-ra,i/roo), 

apooWoi 

k^  =  ^T^exp{T,JT^-Ta,jTf)exp{Ta,rnHT]),  form  #9,  (7) 

Ai 

h  =  t^~TT 

mix 

We  note  that  the  normalized  rate  constants  axe  a  function  of  the  frozen 
temperature,  T/,  and  the  temperature  perturbation,  (p.  A  is  the  Damkohler 
number  based  on  the  hydrogen  branching  rate  and  the  characteristic  con¬ 
vective  rate  of  the  system.  At  ignition,  we  expect  this  ratio  to  be  0(1). 

Examining  the  equation  governing  the  oxygen  radical,  we  insert  pertur¬ 
bations  to  the  frozen  reactant  concentrations  of  the  form  Ybj  =  Yq^j  —  yo, 
and  Ihj  =  Yh^j  ~yH2  >  where  yoi  and  y^^  are  expected  to  be  small  compared 
to  their  respective  frozen  values.  This  results  in  the  following  equation  for 


dYo  r  ^  .r,  ^  ^ 

~  -A[A:iyH(i^2,/  “  yOi) 

—k2{YH2j  —  yH-i)Yo  —  ki2YHO-iYo]- 

The  coefficient  of  Yo  in  the  reaction  mechanism,  k2YH^j,  is  >  1  for  0(1) 
X  and  0(1)  mass  fractions  of  H2  in  the  fuel  stream.  Therefore  the  rate 
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of  reaction  2  will  be  much  larger  than  the  convective-diffusive  rate  for  any 
scale  of  Yo-  As  a  first  approximation,  we  can  neglect  the  convective  and 
diffusive  terms  and  put  the  oxygen  radical  in  chemical  steady  state.  This 
leads  to  the  following  expression  for  Yb: 

Y  =  fciyg(yb;,/  -yoj  /gN 

k2{YHt,f  -  yHa )  +  kuYnOi 

This  approximation  breaks  down  for  large  positive  x  where  Yjj^j  becomes 
0{lfk2).  The  second  order  correction  to  this  approximation  in  the  reac¬ 
tion  mechanism  is  of  0{Yh I k2YH2,f)  which  is  negligible  compared  to  the 
branching  rate. 

In  a  similar  manner,  we  can  write  the  governing  equation  for  the  mass 
fraction  of  Yoh- 

=  -A  lhYH(Yo,.,  -yo,)  +  h(Ytt,.,  -  yH,)Yo 

ax  dx^ 

-hlYn-tj  -  yH2)YoH  +  2knYHYH02  +  kuYnOiYo] 

(10) 

The  coefficient  of  Yb^  in  the  reaction  mechanism,  fcsYffj,/,  is  much  larger 
than  one  for  a;  of  0(1),  so  we  can  again  assume  OH  to  be  in  chemical  steady 
state,  with  wa  =  Wi  +  wj  +  2ujn  +  (^12-  These  approximations  considerably 
reduce  the  complexity  of  the  problem  because  all  occurrences  of  Yb  and 
Yoh  in.  the  reaction  mechanism  can  be  eliminated.  Differential  equations 
need  to  be  solved  only  for  the  remaining  species. 

To  proceed  with  the  analysis,  the  governing  equations  were  put  in  lin¬ 
ear  form  in  a  manner  similar  to  that  of  Sanchez  et  al.  (1994).  Inserting 
expansions  of  the  form: 

Yh  =  eyH,\  +  e^3/H,2,  Yuo^i  =  eyHOz.D  ^^2  =  Yh^J  - 

Yo:,  =  Yo^j  -  ^yo2.i,  <t>  =  A  =  A/  -t-  eAs,  (11) 

and  gathering  terms  of  like  order  in  e,  the  governing  equations  become  at 

0(6): 

^ -Ai{2hj  -  ksj)Yo,jyH,i,  (12) 

ax  dx^ 
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dyH02A  ,  T  d?yH02,\  a  i.  v- 

-  +  Lch — — —  =  -AiksjYo2,fyH,i, 


dx 


dx^ 


dyo2. 


dx 

d(f>i 


-  +  Lea 


d^yo2.i 

dx^ 


+  h,f)Yo2,fyH,i, 


dx'^ 

dyHi.i  ,  d^VHi,! 


^Hx  ~  +  99^./)^02./2/h.1) 


+ 


dx  '  dx^ 

and  at  O(e^): 

-  -Ab(2A:i./  -  k9.f)Yo2,fyH,x 
Ta.x 


—3AikijYo2,fyH,i, 


dx 


dx^ 


(13) 

(14) 

(15) 

(16) 


I  yr  .  1  ^  kujyHO,.! 

,f  \yH,2Yo2.f  +  Yo2,f^<pyH,x  - yH.iyo^.x  - yn^yo^j  / 

—k9,fiyH,2^02,f  —  j/K.iyoj.i)  +  ^11,/  ^1  —  yHOi.iyH.i  , 


(17) 


where  ql  =  qx  +  Q3  +  2^3  and  k^j  is  the  rate  constant  of  reaction  m  evalu¬ 
ated  at  the  frozen  temperature.  The  0(e)  equations  determine  the  ignition 
Damkohler  number,  A/,  while  the  O(e^)  equation  for  the  hydrogen  radical, 
Eq.  (17),  determines  the  behavior  of  the  ignition  point,  Ab-  The  meaning 
of  the  bifurcation  factor,  Ab,  will  be  discussed  in  a  subsequent  section. 

The  first  order  equations  are  equivalent  to  those  which  would  result  from 
applying  the  chemical  steady  state  assumption  of  the  O  and  OH  radicals 
to  reactions  1,  2,  3,  and  9: 


ZH2  d"  O2  — *  2^7  -|-  ‘IH2O  uj  —  tui, 
H  +  O2  +  M  HO2  +  M  w  =  (jJs. 


(18) 


This  mechanism  clearly  shows  the  importance  of  the  crossover  tempera¬ 
ture  defined  by  2ki  =  kg.  For  Too  >  T*,  the  H  production  rate  is  greater 
than  its  termination  rate,  and  a  chain  branching  explosion  is  possible.  For 
Too  <  T*,  chain  branching  explosion  of  H  can  no  longer  occur,  and  ig¬ 
nition  must  proceed  through  a  different  chemical  process.  The  preceding 
reduced  mechanism  is  the  most  fundamental  mechanism  which  can  be  used 
to  determine  the  hydrogen-oxygen  ignition  limits  for  temperatures  above 


the  crossover  temperature. 
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The  first  order  equation  for  the  hydrogen  radical,  Eq.  (12),  is  a  linear 
eigenvalue  problem  which  determines  the  ignition  Damkdhler  number.  The 
initiation  reaction,  RlOb,  has  been  neglected  in  this  equation.  Because  of  its 
slow  rate,  this  reaction  will  only  be  important  on  the  slowly  burning  branch 
of  the  solution.  Since  Eq.  (12)  is  an  eigenvalue  problem,  adding  or  subtract¬ 
ing  a  small  source  term  such  as  initiation  changes  only  the  behavior  of  the 
solution  up  to  the  ignition  point;  the  ignition  Damkdhler  number  remains 
unaffected.  To  confirm  this,  figure  1  compares  the  maximum  value  of  Yh 
determined  from  the  numerical  solution  of  Eq.  (12)  with  various  initiation 
rates  to  the  case  in  which  initiation  has  been  neglected.  The  behavior  shows 
that  neither  the  ignition  Damkdhler  number  nor  the  behavior  of  the  solu¬ 
tion  above  the  ignition  point  is  affected  by  order  of  magnitude  changes  in 
the  initiation  rate.  After  ignition  occurs,  the  sharp  increase  in  the  hydrogen 
radical  concentration  causes  the  rate  of  the  branching  reaction  to  be  orders 
of  magnitude  greater  than  the  initiation  reaction.  We  can  therefore  ignore 
the  effect  of  the  initiation  reaction  in  determining  A/  and  the  behavior  of 
the  solution  at  the  turning  point. 

Another  point  to  note  is  that  the  frozen  fuel  concentration  does  not 
appear  in  the  first  order  equations.  This  results  from  the  steady  state 
approximation  of  the  O  and  OH  radicals.  Because  the  rate  constants  of 
reactions  2  and  3  are  much  greater  than  that  of  reaction  1  these  reactions 
do  not  limit  the  ignition  process.  Since  these  are  the  only  reactions  depen¬ 
dent  on  the  fuel  concentration.  A/  is  independent  of  the  H2  concentration. 
This  is  in  agreement  with  numerical  calculations  which  show  that  the  pres¬ 
sure  and  temperature  ignition  limits  are  insensitive  to  the  H2  concentration 
(Balakrishnan  et  al.,  1995). 

Similarly,  the  higher  order  perturbation,  i/Hj,  was  also  neglected  in  Eq. 
(17).  Depletion  of  H2  in  the  reaction  zone  will  have  a  weak  effect  on  the 
ignition  behavior  because  of  the  weak  dependence  on  the  H2  concentration 
in  the  first  order  equations.  The  effect  of  hydrogen  depletion  in  the  reaction 
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zone  is  captured  only  for  frozen  H2  concentrations  small  enough  to  cause 
departure  from  the  chemical  steady  state  approximations. 

Ignition  Criterion 

In  proceeding  with  the  asymptotic  analysis  to  determine  the  ignition  crite¬ 
rion,  two  possible  limiting  processes  can  be  exploited.  First,  we  can  pursue 
the  large  activation  energy  limit,  which  is  commonly  used  in  problems  with 
one-step  chemistry  (e.g.  Lihan,  1974).  However,  due  to  the  large  diffusivity 
of  the  hydrogen  radical,  it  becomes  difficult  to  match  solutions  between  the 
various  asymptotic  zones  (See  Appendix  A).  To  avoid  this  difficulty,  we 
instead  use  the  Lewis  number  of  the  hydrogen  radical  as  a  small  expansion 
parameter  as  in  Sanchez  et  al.  (1994).  When  the  Lewis  number  is  treated 
as  a  small  parameter,  we  find  that  the  solution  domain  is  divided  into  three 
zones.  In  a  zone  of  width  0(\/Le^)  around  the  stagnation  plane,  the  frozen 
temperature  profile  decays  from  Too  to  T_oo  and  the  frozen  oxygen  profile 
decays  from  one  to  zero.  Upstream  of  the  stagnation  plane  in  the  oxygen 
jet,  the  temperature  and  oxygen  mass  fraction  are  constant  equal  to  Too  and 
1  respectively.  Similarly  in  the  hydrogen  jet,  the  temperature  and  oxygen 
mass  fraction  are  constant  equal  to  T_oo  and  0  respectively.  A  schematic  of 
the  asymptotic  structure  is  shown  in  figure  2.  Ignition  is  expected  to  occur 
in  the  high  temperature  stream  due  to  the  non-zero  oxygen  concentration 
and  the  exponentially  larger  reaction  rates  in  the  high  temperature  stream. 

A  difficulty  in  the  analysis  is  that  equations  (12)-(17)  are  not  valid  for 
a:  1  due  to  the  breakdown  of  the  steady  state  approximations.  To  match 

the  solution  to  the  high  temperature  boundary  through  this  zone,  two  lim¬ 
iting  approaches  were  used.  The  first  approach  assumed  that  the  chem¬ 
ical  steady  state  approximations  remain  valid  to  the  upstream  boundary. 
This  method  overestimates  the  reaction  rates  in  this  zone  and  thus  gives 
a  lower  limit  to  the  ignition  Damkohler  number.  The  second  method  as¬ 
sumed  that  reaction  was  frozen  for  some  x  ^  1.  These  two  methods  give 
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suitable  bounds  to  the  influence  of  reaction  in  the  large  x  zone  on  the  ig¬ 
nition  Damkohler  number.  Appendix  B  compares  the  difference  between 
these  two  approaches.  We  found  that  if  the  steady  state  assumptions  hold 
for  x  1  then  the  difference  between  these  two  approaches  is  negligible. 

The  equation  for  t/h.i  is  uncoupled  from  the  equations  governing  the 
remaining  species,  and  is  solely  responsible  for  determining  the  ignition 
Damkohler  number.  Inserting  expansions  of  the  form  A/  =  /^p-\-y/LeH^P+ 
•  and  t/h.i  =  2/h,i  +  \/^yH,i  +  LeHy^H\  •  •  •  equation  (12) 


ajid  matching,  we  were  able  to  find  an  analytic  expression  for  A/.  To  first 
order,  Yh  simply  diffuses  through  the  inner  zone,  and  as  we  expected,  the 
ignition  criterion  is  determined  by  the  chemistry  in  the  high  temperature 
stream.  The  first  order  expression  for  A/  is  given  by: 


v^  =  -r 


1  —  (2fei^oo  —  fc9^eo)Aj  jY  ^ 


.  (19) 


where  kx^oo  and  are  equal  to  kx  and  kg  evaluated  at  Too,  and  T  is  the 
gamma  function.  Given  that  A:i^oo  is  equal  to  one,  Aj^^  can  be  calculated 
to  be  1.425/(2  -  Ifcg.oo)-  Thus  for  0(1)  values  of  the  effective  branching 
rate,  (2  —  ifcg.co))  the  ignition  Damkohler  number  is  0(1)  as  we  anticipated. 
However,  for  smaller  effective  branching  rates  A/  can  become  1.  As  A/ 
becomes  large,  the  effect  of  reactions  which  initially  had  very  slow  rates  will 
be  amplified  and  cause  the  analysis  to  breakdown.  Indeed,  this  is  how  the 
transition  to  the  third  ignition  limit  occurs.  Peroxyl  path  reactions  become 
important,  and  ignition  no  longer  occurs  simply  through  hydrogen  radical 
branching.  Transition  to  the  third  limit  is  not  captured  by  this  analysis. 

Further  matching  resulted  in  an  integral  expression  for  Aj  '  which  was 
simplified  using  the  activation  energy  as  a  large  expansion  parameter: 

.  (2)  _ ■\/27rAj  ^(2  —  ^9,00) _ 

^  ”  ’F  ([1  -  (2  -  A:9.oo)Ai')]  /2)  -  ^  (l  -  (2  -  ji)  (2  -  ’ 

(20) 

where  ’F  is  the  digamma  function  and  Xr  is  determined  from  the  relation 
exp(-Xr/2)/v^x,.  =  Tl,fTa,x0-  Physically,  Xr  is  the  location  of  the  ex- 
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ponential  decay  of  the  branching  rate  due  to  the  temperature  gradient,  in 
nondimensional  lengths  of  y^A/po CpO.  The  source  of  this  correction  can  be 
attributed  to  the  exponential  decay  of  the  branching  rate  at  y/LeHXr  instead 
of  at  the  stagnation  plane  as  the  first  order  analysis  assumes. 

Combining  the  previous  expressions  results  in: 


A7  = 


1.425  0.97^/L^Xr 

(2  -  k^,oo)  {2-h.coy 


+  0{LeH). 


(21) 


This  expression  has  been  compared  to  the  ignition  Damkohler  number  cal¬ 
culated  from  equation  (12)  numerically  and  found  to  be  accurate  to  within 
3%  over  most  of  the  parameter  range.  We  can  explain  some  fundamental 
properties  of  hydrogen-oxygen  ignition  using  this  expression.  Rewriting  it 
in  terms  of  dimensional  quantities  gives: 


yo 


Wo 


-Ta.l/To, 


pIys. 


-AgT^  =  ac7(  1.425  -h 


0.97\/  LejjXr 
2  ^9,00 


).  (22) 


This  expression  clearly  shows  that  hydrogen-oxygen  ignition  at  tempera¬ 
tures  above  the  crossover  temperature  is  determined  by  a  balance  between 
branching  and  termination  of  the  hydrogen  radical  in  the  high  temperature 
stream  with  a  convective-diffusive  loss  term  which  is  proportional  to  the 
strain  rate. 

Using  this  ignition  criterion,  we  can  determine  the  ignition  pressure  as 
a  function  of  Too  for  a  given  strain  rate,  as  shown  in  figure  3.  These  curves 
behave  in  a  manner  analogous  to  the  “first”  and  “second”  limits  of  the  homo¬ 
geneous  hydrogen-oxygen  system.  In  the  first  limit,  the  pressure  is  low  such 
that  the  termination  rate,  R9,  which  is  a  three  body  reaction,  becomes  slow 
relative  to  the  branching  reaction,  Rl.  In  equation  (22),  we  can  ignore  the 
termination  term  and  the  equation  reduces  to  poOxp(  Ta,i/Too)  =  constant. 
Therefore  in  the  first  limit,  for  a  constant  strain  rate  a  decrease  in  pres¬ 
sure  must  be  accompanied  by  an  increase  in  temperature.  In  the  second 
limit,  the  branching  and  termination  rates  nearly  balance  so  we  can  write 
2Aiexp(-Ta,i/roo)  =  PoA^T^IWaxt,  which  is  simply  the  expression  for  the 
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crossover  temperature.  This  limit  is  insensitive  to  variations  in  the  strain 
rate  due  to  the  small  net  branching  rate,  (2  -  fcg.oo)-  This  can  be  seen  by 
taking  the  derivative  of  the  first  order  expression  for  A/: 

dD  _ 

dAi  ~  1.425’ 

where  D  is  the  net  branching  rate,  (2  -  A:9,oo).  In  the  second  limit,  the  net 
branching  rate  is  small,  therefore  large  changes  in  the  strain  rate  can  be 
balanced  by  small  changes  in  the  net  branching  rate. 

The  ignition  pressures  and  temperatures  predicted  by  this  expression 
show  good  qualitative  agreement  with  the  numerical  calculations  using  full 
chemistry.  However,  the  ignition  temperatures  are  as  much  as  80K  below 
the  numerical  ignition  temperatures  (Kreutz  and  Law,  1995).  The  differ¬ 
ence  could  possibly  be  attributed  to  simplifications  such  as  the  assumption 
of  constant  transport  coefficients,  or  assumptions  regarding  Lti  or  Cp.  In 
any  case,  the  real  value  of  the  asymptotic  approach  is  its  ability  to  reduce 
this  complex  phenomenon  to  an  analytic  expression  which  can  be  used  to 
understand  the  physics  of  the  ignition  criterion. 

Ignition  Behavior 

The  bifurcation  factor,  A^,  determines  how  the  solution  will  behave  as  the 
concentration  of  the  hydrogen  radical,  e,  increases.  When  A^  is  negative, 
increases  in  e  result  in  a  decrease  in  A.  Therefore,  there  is  a  turning  point 
in  the  steady  state  solution  and  the  system  will  ignite  abruptly.  When  Ab 
is  positive,  increases  in  e  result  in  an  increase  in  A.  In  this  case,  the  system 
will  monotonically  transition  to  an  ignited  state.  The  bifurcation  factor  can 
be  solved  for  from  Eq.  (17)  or  determined  using  the  following  compatibility 
condition  which  results  from  integration  of  the  self-adjoint  form  of  that 
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equation: 


Ab  =  ^iVna  2^1./  +  yo2.i  +  yo2j  ) 


-ksjyoi.i  -  ^11./  ( 1  ~ 


ho/ 


J/HOz.l 


(24) 


//" 

J  —CO 


hi.f, 

vlA'^h.f  -  hjU^'^dx. 

) 

Expanding  for  small  Lewis  number  further  simplifies  this  expression.  Due 
to  the  small  diffusivity  of  O2  and  HO2  relative  to  H,  not  only  is  the  oxy¬ 
gen  concentration  exponentially  small  in  the  low  temperature  stream,  but 
the  HO2  radical  concentration  is  also  exponentially  small.  Therefore,  the 
contribution  of  the  integral  from  —00  to  —x  of  0{s/LeH)  can  be  neglected. 
Additionally,  as  a  first  approximation  the  O^s/Lch)  contribution  from  the 
zone  around  the  stagnation  plane  can  also  be  neglected.  Our  final  sim¬ 
plification  is  to  normalize  yuOi,!^  l/Oj.ij  coefficient  of  j/h.i 

in  equations  (13)-(15)  respectively.  This  results  in  the  following  simplified 
expression  for  the  bifurcation  factor: 

Ab  =  J  ^iy]i,\  +  q9h,oo)<f>  +  (1  +  -  h./yoi.i 

^  "  h2,oo 


no, 00 


nijoo 


^11,00  /  ^2^cxi^H2tS 


^s.ooS/bOj.i  ^  ®  ^  dx 


/  r  yH,i[^h,f-hj]/''''dx. 

J —00 


(25) 


Insight  into  the  controlling  mechanism  of  the  ignition  behavior  can  be  gained 
by  examining  the  terms  of  this  equation.  The  coefficient  of  <f>  in  this  equa¬ 
tion  represents  the  effect  of  heat  release  on  the  bifurcation  factor.  Because 
of  the  large  activation  energy  of  the  branching  reaction,  heat  addition  to  the 
system  causes  the  rate  of  reaction  1  to  exponentially  increase  while  the  ter¬ 
mination  reaction  remains  relatively  insensitive  to  heat  addition.  Therefore, 
the  net  branching  rate  increases  causing  the  solution  to  exhibit  a  turning 
point.  The  magnitude  of  this  effect  is  proportional  to  the  heat  release  from 
the  branching  and  termination  steps  and  to  the  activation  temperature.  Its 
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dependence  on  the  activation  energy  reflects  the  nature  of  the  Arrhenius 
term.  Large  activation  energy  relative  to  the  flow  temperature  indicates 
that  only  a  small  fraction  of  the  possible  branching  collisions  have  enough 
energy  to  react.  In  this  case,  the  branching  rate  will  be  very  sensitive  to 
heat  release.  At  higher  system  temperatures  a  greater  percentage  of  the 
collisions  react  and  the  system  is  less  sensitive  to  heat  release. 

The  second  term  in  the  compatibility  condition  is  due  to  the  effect  of 
oxygen  depletion  on  the  ignition  process.  Depletion  of  oxygen  causes  a  de¬ 
crease  in  the  net  branching  rate.  The  magnitude  of  this  effect  is  proportional 
to  the  amount  of  oxygen  depletion,  (fci  -h  Aig)  as  found  from  Eq.  (14),  times 
the  net  branching  rate,  (2fci  -  fcg).  The  sign  of  this  bifurcation  parameter  is 
positive,  meaning  that  oxygen  depletion  will  reduce  the  net  branching  rate 
and  lead  to  a  monotonic  ignition  phenomenon. 

The  remaining  terms  in  the  bifurcation  factor  are  due  to  the  effect  of 
HO2  chemistry.  These  terms  are  all  proportional  to  kg  because  this  reac¬ 
tion  determines  the  production  rate  of  HOg.  Because  reaction  9  is  a  three 
body  reaction,  its  production  rate  of  HO2  is  slow  at  low  pressures.  There¬ 
fore,  reactions  involving  HO2  are  unimportant  at  low  pressures.  At  higher 
pressures,  however,  production  of  HO2  is  not  negligible,  and  HO2  chem¬ 
istry  can  have  several  significant  effects.  Reaction  11,  H  +  HO2  T-OH, 
is  an  additional  branching  step  for  the  hydrogen  radical,  while  reaction  10, 
H  +  HO2  H2  +  O2,  is  a.  termination  step.  The  ratio  of  the  rates  of  these 
reactions  is  nearly  independent  of  temperature  with  fcio/^ii  being  approx¬ 
imately  1/6.  Lumping  these  two  reactions  together,  the  net  effect  is  the 
production  of  an  additional  5/6  of  a  hydrogen  radical  for  every  reaction  be¬ 
tween  H  and  H02-  Thus  as  HO2  is  produced,  there  is  an  additional  source 
of  branching  for  the  hydrogen  radical.  This  accelerates  the  ignition  process 
and  leads  to  a  turning  point  in  the  steady  state  solution. 

The  final  term  is  due  to  the  effect  of  reaction  12  on  the  steady  state 
approximation  for  the  oxygen  radical.  Reaction  12  competes  with  reaction 
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2  for  oxygen  radicals: 


R2  O  +  H2  — ^  OH  +  -ffj 
R12  0  +  H02^0H  +  02^ 


(26) 


While  reaction  2  is  a  branching  step,  reaction  12  is  only  a  propagation  step. 
At  low  HO2  concentrations  most  of  the  radicals  are  consumed  by  reaction 
2.  However  at  larger  HO2  concentrations  more  0  radicals  are  consumed 
by  reaction  12,  eind  the  net  branching  rate  decreases.  The  bifurcation  pa¬ 
rameter  is  a  ratio  of  the  rates  of  these  two  pathways.  The  sign  of  this 
term  indicates  that  as  H O2  accumulates  there  will  be  a  reduction  in  the  net 
branching  rate.  The  dependence  of  reaction  2  on  the  H2  concentration  is 
also  apparent.  Increasing  the  hydrogen  concentration  increases  the  rate  of 
reaction  2  and  thus  diminishes  the  importance  of  reaction  12  in  determining 
the  bifurcation  factor. 

At  different  pressures,  temperatures,  and  fuel  concentrations  any  one 
of  these  nonlinear  effects  can  control  the  ignition  process.  In  figure  4,  the 
bifurcation  parameters  associated  with  heat  release,  oxygen  depletion,  and 
reactions  10  and  11  calculated  along  the  ignition  curve  of  a  =  1000s“^  are 
plotted  as  a  function  of  temperature.  The  effect  of  reaction  12  will  be  ex¬ 
amined  subsequently  because  of  its  dependence  on  the  fuel  concentration. 
In  the  second  limit,  the  termination  rate  is  0(1),  and  the  ignition  mech¬ 
anism  is  dominated  by  the  effect  of  the  reaction  H  -{■  HO2  20 H.  In 
the  first  ignition  limit,  kg  becomes  <C  1,  and  the  effect  of  heat  release  be¬ 
comes  the  dominant  effect.  At  very  low  pressures,  the  ignition  temperature 
increases,  and,  as  discussed,  the  effect  of  heat  release  diminishes.  Oxygen 
depletion  then  competes  with  heat  release  in  determining  the  ignition  be¬ 
havior.  Figure  4  demonstrates  the  various  responses  that  can  be  exhibited 
by  the  hydrogen-oxygen  system  at  different  ignition  pressures  and  temper¬ 
atures. 

The  bifurcation  parameter  associated  with  reaction  12  and  the  param¬ 
eter  due  to  reactions  10  and  11  are  plotted  along  the  ignition  curve  of 
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a  =  10005“^  for  various  fuel  concentrations  in  figure  5.  This  figure  shows 
that  these  two  bifurcation  parameters  exhibit  similar  temperature  depen¬ 
dence.  In  addition,  for  H2  concentrations  in  the  reaction  zone  for  which  the 
steady  state  assumption  still  holds,  i.e.  1,  the  effect  of  reaction 

12  is  always  dominated  by  reaction  11.  Only  for  small  fuel  concentrations 
will  reaction  12  have  a  significant  effect. 

Figure  6  shows  the  bifurcation  factor  calculated  from  the  compatibility 
condition  along  the  ignition  curve  of  a  =  1000s  for  an  incoming  stream 
of  pure  hydrogen.  Only  the  bifurcation  factor  calculated  along  the  first 
limit  is  shown.  In  the  second  limit,  the  bifurcation  factor  is  negative  with 
a  magnitude  >  1  due  to  the  strong  effect  of  reaction  11.  This  figure  shows 
that  around  2000K  the  bifurcation  factor  crosses  zero.  This  is  in  agreement 
with  figure  4  which  shows  that  at  this  temperature  the  effects  of  heat  re¬ 
lease  and  oxygen  depletion  are  of  the  same  order.  Above  this  temperature, 
the  bifurcation  factor  becomes  positive,  and  the  system  no  longer  ignites 
abruptly. 

Finally,  the  ignition  temperatures  and  pressures  are  plotted  in  figure  7 
with  curves  delineating  the  transition  between  ignition  through  HO2  branch¬ 
ing,  ignition  through  heat  release,  and  the  monotonic  transition  caused  by 
oxygen  depletion  in  the  reaction  zone.  These  curves  were  calculated  by 
equating  the  relevant  bifurcation  parameters.  In  each  of  the  regimes,  the 
ignition  criterion  is  still  determined  by  a  balance  between  branching,  termi¬ 
nation,  and  convective-diffusive  loss  of  the  hydrogen  radical,  but  the  mecha¬ 
nism  by  which  the  system  ignites  is  different.  In  regime  I,  ff-b  JJ02  — >  20 H 
accelerates  the  net  branching  rate  and  causes  the  system  to  ignite  abruptly. 
In  regime  II,  heat  relea.se  accelerates  the  branching  rate  by  its  effect  on 
the  Arrhenius  term  of  reaction  1,  and  again  the  system  abruptly  ignites. 
Regime  III  is  the  monotonic  transition  region.  As  the  Damkohler  number  is 
increased,  diffusion  of  oxygen  limits  the  system,  and  the  system  will  mono- 
tonically  transition  to  an  ignited  state.  In  these  ignition  processes,  once  the 
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system  becomes  unstable  and  begins  to  transition  to  an  ignited  state,  we 
note  that  other  effects  become  important.  As  the  system  ignites  in  regime 
I,  eventually  heat  release  will  also  contribute  to  the  transition  to  a  diffu¬ 
sion  flame.  It  is  only  the  mechanism  by  which  ignition  is  initiated  which  is 
considered  here. 

Conclusions 

Ignition  of  a  hydrogen-oxygen  counterflow  system  for  temperatures  above 
the  crossover  temperature  was  studied  asymptotically  using  eight  elemen¬ 
tary  steps  to  describe  the  high  temperature  ignition  reaction  mechanism. 
For  0(1)  fuel  concentrations,  chemical  steady  state  of  0  and  OH  were  shown 
to  be  justifiable  approximations,  and  a  simplified  reduced  mechanism  was 
derived.  Using  this  reduced  mechanism,  we  were  able  to  verify  the  fact  that 
hydrogen-oxygen  ignition  is  determined  by  a  balance  between  branching, 
termination,  and  transport  of  the  hydrogen  radical.  In  addition,  an  alge¬ 
braic  ignition  criterion  was  found  which  predicted  ignition  limits  analogous 
to  the  first  and  second  limit  of  the  homogeneous  hydrogen-oxygen  system. 
This  criterion  clarified  the  physics  underlying  the  behavior  of  these  limits, 
and  captured  the  dependence  of  the  ignition  limits  on  all  system  parameters. 

Performing  a  bifurcation  analysis  similar  to  that  of  Sanchez  et  al.  (1994), 
the  ignition  behavior  of  the  hydrogen-oxygen  system  was  also  resolved.  We 
found  that  for  temperatures  above  the  crossover  temperature,  abrupt  ig¬ 
nition  can  occur  either  through  interactions  between  H  and  HO2  via  the 
reaction  H  -j-  HO2  20H,  or  through  the  effect  of  heat  release  on  the 
branching  rate.  At  temperatures  much  greater  than  the  crossover  tempera¬ 
ture,  oxygen  depletion  in  the  ignition  zone  can  cause  the  system  to  exhibit 
a  monotonic  transition  to  the  ignited  state. 
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Appendix  A 


When  using  the  activation  energy  as  a  large  expansion  parameter,  the  phys¬ 
ical  domain  can  be  divided  into  three  zones.  In  the  high  temperature 
stream,  there  is  a  zone  in  which  T  =  Too  3-nd  Yoj  =  1.  This  is  followed 
by  a  thin  zone  in  which  the  branching  rate  exponentially  decays  and  then 
a  convective-diffusive-termination  zone.  In  the  upstream  zone,  the  solution 
for  the  hydrogen  radical  is  the  same  as  in  the  small  Lewis  number  case, 
g-LeHxV4(7(x/2  -  AjX^Len),  where  x  is  now  made  nondimensional  by  the 
thermal  length,  and  U  is  the  parabolic  cylinder  function  (Abro- 

mowitz  and  Stegun,  1965).  In  the  inner  zone,  the  solution  satisfies  the 
equation: 


(  1 

XLen  di  Lea  df 


(Al) 


where  i  and  Xr  is  the  large  expansion  parameter  which  is 

determined  by  the  location  of  the  decay  zone.  The  matching  zone  for  the 
high  temperature  stream  and  the  inner  solution  is  x~'^  «C  (s  —  Xr)  <  ^r- 
To  match  these  solutions,  an  expansion  for  e''^®"®^/'‘C/(l/2  —  A,x^/L^)  is 
needed  around  the  point  Normally,  Xr  is  a  large  number  and  a  large- 
X  expansion  can  be  used,  but  because  the  coupling  Xts/Lch  is  actually 
less  than  one,  this  method  would  be  inappropriate.  Furthermore,  once  the 
matching  condition  was  derived,  equation  (Al)  would  still  need  to  be  solved 
numerically  to  determine  the  ignition  Damkohler  number. 

The  small  Lewis  number  method  (Sanchez  et  al,  1994)  allows  us  to 
avoid  these  difficulties  and  obtain  an  analytic  expression  for  the  ignition 
Damkohler  number.  It  is  also  a  more  appropriate  expansion  parameter 
because  the  large  activation  energy  expansion  parameter  Xr  is  only  slightly 
larger  than  one,  whereas  the  Lewis  number  is  0(10“^). 


22 


Appendix  B 

The  solution  in  the  high  temperature  chemical  steady  state  region  is; 

cie-"'/^C/(l/2  -  A,  a:)  +  C2e-"'/'‘K(l/2  -  A,  x),  (Bl) 


where  U  and  V  are  the  parabolic  cylinder  functions  (Abromowitz  and  Ste- 
gun,  1965).  Assuming  that  this  solution  is  valid  through  the  region  where 
the  steady  state  approximations  degenerate,  we  can  find  a  relation  between 
Cl  and  C2  by  enforcing  the  condition  that  Yh  —  0  aX  the  high  temperature 
boundary.  Using  the  large-x  expansions  for  these  functions  the  following 
equation  is  arrived  at: 

+  C2X^^  =  0,  or  C2/C1  =  (B2) 


where  xj,  is  the  boundary  location.  In  the  limit  of  X5  00,  we  see  that  C2  is 
negligible  and  the  solution  in  the  high  temperature  stream  is  Cie“®^/^U(l/2— 
A,x). 

When  reaction  is  assumed  frozen  in  the  region  where  the  chemical  steady 
state  approximations  degenerate,  we  need  to  perform  matching  between  the 
frozen  region  and  the  chemical  steady  state  region.  In  the  frozen  region  we 
have  the  solution  di  +  (i2crfc(x/  v^),  and  in  the  chemical  steady  state  region 
the  solution  is  again  as  in  equation  (Bl).  Enforcing  the  high  temperature 
boundary  condition  to  the  frozen  solution  in  the  limit  of  the  boundary 
location  going  to  infinity  results  in  the  solution  d2erfc(x/\/2)  in  the  frozen 
zone.  If  we  match  slopes  and  values  between  the  frozen  solution  and  the 
steady  state  zone  solution  at  a  location  where  x  >  1  but  the  steady  state 
approximations  are  still  valid,  we  get: 


c.e-5'=  +  c.*7'^,  (B3) 

-(A-l)xf^)-c,Ax;^-‘.  (B4) 
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Solving  for  the  ratio  C2/C1  gives: 

^  -  l)xf (B5) 

Cl  v27r 

where  x/  is  the  matching  location.  This  ratio  is  again  exponentially  small 
for  Xf  ^  1.  The  solution  in  the  steady  state  region  is  the  same  for  both 
cases,  which  shows  that  the  non-steady  state  region  has  little  importance  in 
determining  the  ignition  criterion.  However ,  if  the  incoming  fuel  concentra¬ 
tion  becomes  <  1,  then  the  steady  state  approximations  will  break  down 
for  X  of  0(1)  and  there  will  indeed  be  an  effect  on  the  ignition  criterion  due 
to  the  breakdown  in  the  steady  state  approximations. 
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Figure  Captions 


Figure  1  Comparison  of  numerical  solutions  for  maximum  Yh  with  and 
without  initiation  reaction. 

Figure  2  Schematic  of  the  frozen  temperature,  oxygen,  and  hydrogen 
profiles  in  each  asymptotic  zone. 

Figure  3  Pressure  and  temperature  ignition  limits  for  various  strain 
rates. 

Figure  4  Bifurcation  parameters  of  reactions  10  and  11,  heat  release, 
and  oxygen  depletion  along  the  first  and  second  ignition  limits 
of  a  =  lOOOs”^. 

Figure  5  Bifurcation  parameters  of  reactions  10  and  11,  and  reaction 
12  along  the  first  and  second  ignition  limits  of  a  =  lOOOs"^ 
for  various  fuel  concentrations. 

Figure  6  The  bifurcation  factor,  Ab,  along  the  first  ignition  limit  of 
a  =  lOOOs”^. 

Figure  7  High  temperature  ignition  regimes  of  the  hydrogen-oxygen 
system. 
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Table  1  Starting  reaction  mechanism  for  the  hydrogen-oxygen  system. 
Rate  constants  are  in  the  form  k  =  AT^  Eal R°T).  Units  are  moles, 

cm^,  seconds,  Kelvin  and  kcal/mole. 


no. 

reaction 

A 

n 

E. 

R1 

H  +  O2  0  +  OH 

1.92E-t-14 

0.00 

16.44 

R2 

0  +  H2^  H  +  OH 

5.08E-1-04 

2.67 

6.29 

R3 

OH  +  H2^  H  +  H2O 

2.16E+08 

1.51 

3.43 

R9 

H  +  O2  +  M  ^  HO2  +  M 

6.17E-t-19 

-1.42 

0.00 

RIO 

HO2  -\-  H  — >  H2  +  O2 

6.63E+13 

0.00 

2.13 

RlOb 

+  O2  ->  HO2  +  H 

1.93E+14 

0.00 

59.61 

Rll 

HO2  aH  -^OH  +  OH 

1.69E-fl4 

0.00 

0.87 

R12 

HO2  AO  OH  AO2 

l,81E“hl3 

0.00 

-0.40 
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